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ABSTRACT 

Application of multi-resolution analysis adopting wavelets in order to investigate mil- 
lisecond aperiodic X-ray variability of Cyg X-l is presented. This relatively new ap- 
proach in a time-series analysis allows us to study both the significance of any flux 
oscillations localized in time-frequency space, as well as the stability of its period. Us- 
ing the observations of the RXTE/PGA we analyze archival data from Cyg X-l during 
its failed state transition. The power spectrum of such a state is strongly dominated by 
a single Lorentzian peak corresponding to a damped oscillator. Our wavelet analysis 
presented in this paper shows the existence of short lasting oscillations without an 
obvious trend in the frequency domain. Based on this result, we suggest the interpre- 
tation of the dominant components of Cyg X-l power spectrum in all spectral states: 
(i) a power law component dominating a soft state is due to accretion rate fluctua- 
tions related to MRI instability in the cold disk, (ii) the low frequency Lorentzian is 
related to the same fluctuations but in the inner ion torus, and (in) the high frequency 
Lorentzian (dominating the power spectrum during the failed state transition) is due 
to the dynamical pulsations of the inner ion torus. 

The usage of the wavelet analysis as a potential and attractive tool in order to 
spot the semi-direct evidences of accretion process onto black-holes is proposed. 

Key words: accretion, accretion disks - binaries: general - stars: individual: Cyg X-l 
- X-rays: observations - X-rays: stars. 



1 INTRODUCTION 

Cyg X-l (Bowyer et al. 1965) belongs to the class of Galac- 
tic accreting black-hole binaries characterized by a highly 
variable X-ray emission in time-scales covering up to ten or- 
ders of magnitude (e.g. Nowak et al. 1999; Reig, Papadakis 
& Kylafis 2002; Psaltis 2004). The optical companion, HDE 
226868 star (Bolton 1972; Webster & Murdin 1972) was clas- 
sified as 09.7Iab supergiant (Walborn 1973) with the orbital 
period of 5.6 d. The main source of the matter accreting onto 
black hole is the focused stellar wind taking its origin from 
HDE 226868. Mass ranges of the components are 10-32 M 
and 16-55 Mq for the black-hole and the supergiant, respec- 
tively (Gies & Bolton 1986; Herrero et al. 1995; Ziolkowski 
2005). 

The power spectrum density (PSD) of the X-ray emis- 
sion is broad-band and aperiodic. It was frequently char- 
acterized as a broken power-law (Sutherland, Weiskopf & 
Kahn 1978; Belloni & Hasinger 1990; Nowak et al. 1999; 
Revnivtsev, Gilfanov & Ctmrazov 2000), with occasional 



E-mail: paulo@camk.edu.pl 



narrow quasi-periodic features (QPO) (Angelini, White & 
Stella 1992). Higher quality data allowed for more sophis- 
ticated parameterization of the PSD through a number of 
broad Lorentzians and a power-law component (e.g. Nowak 
2000; Pottschmidt et al. 2003, hereafter P03). The lightcurve 
is not stationary. In short time-scales the variance is pro- 
portional to the mean luminosity of the source (Uttley & 
McHardy 2001; Gleissner et al. 2004) so the normalized 
power spectrum is broadly used to characterize the source 
(Miyamoto et al. 1992 and most of the subsequent papers). 
In long time-scales also the shape of the PSD changes as it 
depends significantly on the spectral state of the source (e.g. 
Cui et al. 1997a; 1997b). 

Although frequently used, Fourier-based technique is 
actually not suitable for the analysis of the aperiodic vari- 
ability. The Fourier analysis uses a fixed sinusoidal basis 
functions to decompose a signal. This method therefore pro- 
vides excellent frequency resolution of persistent features 
but it fails when the signal is highly aperiodic and existing 
modulations are needed to be localized in time-frequency 
space. 

To overcome this problem one needs to use a multi- 
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Figure 1. The lightcurve of Cyg X-l as observed by the 
RXTE/PCA on 1999.12.05, showing a 60 s sequence of data se- 
lected by us for deeper studies using wavelet analysis. 



resolution approach, e.g. by wavelet analysis (WA). Because 
the wavelet decomposition gives an idea of both the local 
frequency content of a time-series and the temporal dis- 
tribution of these frequencies, its application in astronom- 
ical signal processing is increasing rapidly (e.g. Szatmary, 
Vinko & Gal 1994; Escalera & MacGillivray 1995; Frick et 
al. 1997; Aschwanden et al. 1998; Barreiro & Hobson 2001; 
Freeman et al. 2002; Irastorza et al. 2003). Interestingly, 
so far the application of wavelets in order to study vari- 
ability of X-ray binary systems and active galactic nuclei 
was marginal. Scargle et al. (1993) used them to exam- 
ine QPO and very low- frequency noise for Sco X-l whereas 
Steiman-Cameron et al. (1997) supplemented Fourier detec- 
tion of quasi-periodic oscillations in the optical lightcurves 
of GX 339-4 by scalegrams, a wavelet technique. Fritz & 
Bruch (1998) applied this technique to the analysis of the 
optical luminosity variations in cataclysmic variables while 
Liszka, Pacholczyk & Stoeger (2000a) analyzed the ROSAT 
lightcurve of NGC 5548 in short time-scales. 

In this paper we select a very special observation of 
Cyg X-l performed on 1999.12.05. It was examined by P03 
and the state of the source was described as failed state tran- 
sition. The power spectrum at that time was particularly 
simple since the dominant part of the power was contained 
in a single Lorentzian peak. The aim of our paper is to verify 
with the wavelet analysis whether indeed a simple damped 
oscillator well represents the time-resolved lightcurve prop- 
erties. 

Our paper is organized as follows. Section 2 provides in- 
formation on our data selection and reduction. In Section 3 
we describe most essential aspects of the time-frequency 
analysis, starting from a brief discussion of computing X- 
ray power spectra and associated time-resolution problem. 
Next, we provide a description of calculation of a continu- 
ous wavelet transform in the frame of our interest. Section 4 
contains the analysis of Cyg X-l lightcurve whereas a discus- 
sion of results, including the lightcurve simulations, is the 
subject of Section 5. We conclude in Section 6. 



2 DATA SELECTION AND REDUCTION 

To probe the object variability in time-scales of seconds we 
use the observation by Proportional Counter Array on board 
of the Rossi X-ray Timing Explorer (RXTE/PCA; Bradt, 
Rothschild & Swank 1993) taken from the public archive of 
the HEASARC. We select the PCA data set of 40099-01-24- 
01 (1999.12.05) on which we concentrate our investigation. 
The data were reduced with the LHEASOFT package v. 
5.3.1 (May 2004) applying a standard data selection, i.e. 
the Earth elevation angle > 10°, pointing offset < 0°.01, the 
time since the peak of the last South Atlantic Anomaly SAA 
> 30 min and the electron contamination < 0.1. The number 
of Proportional Counting Units (PCUs) opened during the 
observation was equal 3. 

We extract the lightcurve in 2.03-13.1 keV energy band 
(channels 0-30) and rebin it with one common bin size 
value of At = 0.025 s. The average count-rate in this set 
is x = 6061.8 cts s _1 and the fractional rms 21.6 %. For 
needs of computation of the Fourier power spectra, we use 
POWSPEC software with an applied Poisson noise level sub- 
traction. For purposes of wavelet analysis, the IDL software 
provided by Torrence & Compo (1998) (hereafter TC98; see 
Acknowledgments for the URL address) was run. 



3 TIME-FREQUENCY ANALYSIS 

Computation of the Fourier power spectral density is an 
ideal tool for detecting periodic fluctuations in the time- 
series. Therefore, in the absence of major problems caused 
by a window pattern (aliasing) and at high signal-to-noise 
ratio, a strong peak in the Fourier spectrum has an imme- 
diate interpretation (e.g. Gray 1976). 

However, the Fourier analysis fails when time evolution 
of spectral features has to be taken into account. In order 
to capture this effect one needs to apply a time-frequency 
representation of the signal. The simplest solution to this 
problem is the application of Windowed Fourier Transform 
(see e.g. Cohen 1995; Flandrin 1999 for review). More ad- 
vanced methods include inter alia Wigner-Ville distribution, 
the reassignment method or the Gabor spectrogram (see Vio 
& Wamsteker 2002 and references therein for a nice review) . 
In the present paper we apply the alternative approach: the 
wavelet analysis. 



3.1 Fourier power spectrum 

For a reference, we calculate the PSD of the studied 
lightcurves using a standard approach (see van der Klis 
1989 for review). We use the normalization of Miyamoto 
et al. (1992) which gives the periodogram in units of 
(rms/mean) 2 /Hz. Therefore, the integrated periodogram 
yields the square of fractional rms of a lightcurve. We sub- 
tract the noise as in Vaughan et al. (2003), omitting the 
correction for dead-times (see Zhang et al. 1995 for the dis- 
cussion of the dead-time effect), as our data sets are rather 
short (20 s) and thus the uncertainty in the determination 
of the PSD dominates. 
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3.2 Wavelet analysis 

A basic concept standing for its application to the time- 
series was to analyze the signal at different frequencies with 
different resolutions (see e.g. Farge 1992; Daubechies 1992 
and TC98 for good summary on wavelets). Contrary to the 
Fourier transform, the wavelet analysis makes use of a set of 
functions, wavelets, which are localized in scale-time space 1 . 
Several standard shapes are broadly used, depending on the 
subject under study. 



3.2.1 Wavelets 

The basic concept is a mother wavelet function, ip(t), which 
depends on time-parameter t. A wavelet's family is therefore 
generated by scaling and shifting of ip{t) as follows: 



1> 



, b G 5ft, a > 



(1) 



where a and b are factors standing for scaling and time- 
shifting, respectively. 

There are three most commonly used wavelets in the 
astronomical time-series analysis: (a) a Marr or Mexican 
Hat: 



V>Marr(t) = (1 — t 2 )e ' /2 

(b) a Morlet wavelet: 

V , Morlot(t) = 7T 



1/4 i2nu t^ — t A /2 



(2) 
(3) 



which has the complex sinusoidal waveform confined by the 
Gaussian bell envelope (Grossman & Morlet 1984), and (c) 
a Paul wavelet: 



V'Paul(i) = 



Vvr(2m)! 



it) 



-(m+l) 



(4) 



where m denotes an order (Paul 1984; Combes, Grossmann 
& Tchamitchian 1989). These functions are presented in 
Fig. [5] Application of complex wavelets instead of real- 
valued ones is required to capture information on the ampli- 
tude and phase simultaneously. On the other hand, a real- 
valued functions are useful for isolation of positive and neg- 
ative modulations as the separate peaks in the time-scale 
plane, making a clear discrimination of the sharp features 
and the signal discontinuities possible to resolve. 

In our research we are mainly interested in detection 
and studies of aperiodic flux structures, therefore the Morlet 
wavelet (with 2-kvo — 6) will be used throughout the paper. 
As it oscillates due to a term of e lt and owns a complex form 
it becomes a desired probing tool in tracing and quantifying 
quasi-periodic modulations at different scales. 
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Figure 2. Examples of mother wavelets: (a) Marr wavelet; (b) 
Morlet wavelet; (c) Paul wavelet. The real and imaginary parts 
of complex wavelets were plotted by solid and dashed lines, re- 
spectively. 



3.2.2 Wavelet power spectrum 

For discretely defined time-series Xk, k = 0, iV bs — 1, the 
wavelet transform can be denoted in its discrete form as: 



k'=0 



X k >4> 



(k' - k)At 



(5) 



In the wavelet analysis a term scale is used instead of frequency 
which is reserved for the Fourier analysis. 



where At is the sampling time of a lightcurve, k represents a 
localized time index and a factor of \J At/a ensures that the 
wavelet contains the same energy everywhere in time-scale 
space of Wk{a) coefficients. In practice, there is no need to 
use formula \Q for the calculation of wavelet transform. It 
can be expressed in terms of the inverse Fourier transform 
of the product of the Fourier transforms of the signal and 
wavelet function as given by equation lAlU and therefore 
allow to speed up the calculations. 
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Frequency, v (Hz) 

Figure 3. Fourier power spectrum of Cyg X-l on Dec 5, 
1999. Solid line represents a broad Lorentzian peak fitted by 
Pottschmidt et al. (2003). 



A wavelet power spectrum can be simply denned as nor- 
malized square of the modulus of the wavelet transform: 



W = (,\w a ,b\ = £\wk(a m )\ 



(6) 



A choice of a normalization factor can be done arbitrarily 
where the most convenient way is to set £ = cr -2 , i.e. equal 
to the inverse square of the lightcurve variance. Here, a 2 
represents the expectation value of the wavelet transform for 
a white-noise process at each scale a m and a time location k 
(TC98). The local wavelet power spectrum W with £ = a~ 2 
is distributed as xi> i- e - as X 2 distribution with two degrees 
of freedom, therefore it allows to determine the contours at 
given confidence level for every flux oscillation. A reality 
of any peak in the wavelet map is tested against certain 
background spectrum and thus any peak can be accepted or 
rejected at earlier assumed significance level. 

Due to a finite duration of a lightcurve some artifacts 
may appear at the edges of the wavelet maps. To force them 
to be negligible one may introduce a cone of influence de- 
fined as a region of the wavelet map where a wavelet power 
in the vicinity of signal discontinuity decreases by the factor 
of e- 2 (TC98). 

2-D information about the signal in the wavelet map 
also allows us to compare it to the results of Fourier power 
spectrum. It can be achieved by performing an integration 
of the wavelet spectrum over time which is known under 
the name of a global wavelet (power) spectrum or a scale- 
gram technique (Scargle et al. 1993). Its discrete form can 
be denoted as: 

k 2 



G{fln 



N 



k — ki 



where < ki < &2 < iV bs — 1 and £ denotes the normaliza- 
tion factor. Here, k\ and hi ought to be chosen to lie outside 
the cone of influence for particular scale a m . Since the nor- 
malization can be arbitrary, one can choose it the same as 
for the wavelet power spectrum (e.g. £ = °" -2 ) an d there- 
fore determine the proper significance levels (see TC98 for 
details). If so, when a comparison of J2j with the calculated 
Fourier power spectrum is required, it is possible to adopt 
a normalization factor of the wavelet spectrum £ = 2At/x 2 . 
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Figure 4. Exemplary power spectra calculated for the single 
consecutive data streams of duration 20 s corresponding to the 
PCA lightcurve intervals: (a) 1090-1110 s, (b) 1110-1130 s, (c) 
1130-1150 s. For purposes of a comparison with the global wavelet 
power spectra in Fig. [F] the subtraction of a noise level was not 
applied. 



Thus, the global wavelet power spectrum will have the units 
of (rms/mean) 2 /Hz. 

A relation between the wavelet scale a and the Fourier 
frequency v depends on our choice of the wavelet func- 
tion. For the Morlet function {2-nuo = 6) it holds i^Morict — 
(1.03a) -1 , for the Marr wavelet KMarr — (3.97a) -1 whereas 
for the Paul wavelet (m = 6) ^p au i ~ (0.97a) -1 . Hereafter, 
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y-axes for all of our wavelet maps will be expressed as fre- 
quency, given by the above transformations. 



4 RESULTS 

The selected data set covers 1364 s. In Fig. 3 we show the 
corresponding PSD, together with the dominant Lorentzian 
component fitted by P03. The shape of the Lorentzian com- 
ponent was described as: 



L 2 {u) 



2R\Q 2 v 2 



(8) 



and the best fit values of the parameters were: amplitude 
R 2 = 0.188 ± 0.004, center frequency v 2 = 1.769±o.oto Hz > 
and the quality factor Q 2 = 0.395lg'Q3g. This power spec- 
trum component has a maximum at the frequency 2.856 Hz, 
and a characteristic decay time-scale, r 2 = Q 2 /(iru 2 ), is 
given by t 2 = 0.07 s. 

The presented PSD shows relatively small errors since 
the whole time sequence was used. The whole lightcurve was 
divided into intervals of 25.6 s, the PSD was calculated for 
each part separately, and finally averaged PSD was plotted, 
as customary done for such data. 

When single PSD spectra are not averaged, they show 
a considerable scatter. This scatter is expected from the \ 2 
distribution of the individual PSD values. For the purpose of 
further discussion of the wavelet maps, we show (see Fig. 4) 
the examples of the PSD obtained for the single fragments 
of the lightcurve of duration 20 s. Errors are large but the 
Lorentzian shape is still visible. The selection of 60 s long 
sequence from whole duration of our PCA time-series was 
done by us randomly, here corresponding to the interval be- 
tween 1090-1150 s (Fig.|TJ. 



Table 1. Properties of the localized peaks in the wavelet maps 
for the assumed significance level of 90%. 

Sequence Cyg X-l Simulated lightcurve (equation |5J 
Number of localized peaks 



Seq. 1 


IS 


11 


Seq. 2 


16 


16 


Seq. 3 


17 


12 




Activity level, k, in % 


Seq. 1 


63.6 


51.2 


Seq. 2 


59.9 


59.5 


Seq. 3 


62.3 


67.9 




Median of peak duration, 


AT P , in [s] 


Seq. 1 


0.54 


0.54 


Seq. 2 


0.54 


0.45 


Seq. 3 


0.27 


0.54 




Median of peak width, 


A log Up 


Seq. 1 


0.12 


0.13 


Seq. 2 


0.14 


0.11 


Seq. 3 


0.14 


0.15 




Median of peak frequency, 


Up, in [Hz] 


Seq. 1 


3.7 


3.9 


Seq. 1 


4.3 


3.8 


Seq. 3 


5.2 


2.8 




Peak width, A log u^ , of the highest peak 


Seq. 1 


0.25 


0.23 


Seq. 2 


0.21 


0.35 


Seq. 3 


0.39 


0.47 


Peak duration, AT^ , of the highest peak in [s] 


Seq. 1 


1.1 


1.2 


Seq. 2 


2.8 


1.6 


Seq. 3 


2.7 


2.7 


Peak frequency, u^ , of the highest peak in [Hz 


Seq. 1 


3.4 


2.2 


Seq. 2 


2.1 


1.7 


Seq. 3 


1.1 


1.5 



4.1 Wavelet maps 

We construct the wavelet maps for the same parts of the 
lightcurve as used in the PSD analysis. The selected band 
of scales corresponds to the Fourier frequency range between 
1 and 10 Hz, where the Lorentzian has the maximum. 

We use the Morlet wavelet, with the standard assump- 
tion 2-7r^o = 6 (see Section ^. 2. lfl . This means that we probe 
the signal with a damped wave performing ~ 3 oscillations. 

The result is shown in Fig.|3 for three consecutive parts 
of the lightcurve. Integrated wavelet spectra are shown to 
the right of each of the sequences. 

In all maps many localized strong peaks are seen, stand- 
ing clearly out of the background. This is even better seen 
in a 3-D plot (see Fig. HJ upper panel). Those peaks domi- 
nate the integrated spectrum although they are present only 
occasionally, mostly between 2-5 Hz. Many lower but still 
significant peaks are also present, so active phase with sig- 
nificant peaks covers more than half of the observed time. 
No single frequency seems to be favored and we do not see 
any particular evolutionary trend. The largest peaks extend 
in time typically for ~ 1 s which means that they are prac- 
tically unresolved. Those which extend for a few seconds 
can be elongated either in the direction of lower or higher 
frequencies, and the shapes are rather irregular. 

We analyzed the distribution and the properties of the 



significant peaks present in the observational data. The re- 
sults are given in Table For the analysis, we took only 
those peaks which are well determined, i.e. detected with 
90% confidence level and not extending beyond the studied 
frequency range. We give there the number of such localized 
peaks in each of the three sequences. We define the activ- 
ity level, k, in each sequence as the fraction of time when at 
least one peak at any frequency is present. We determine the 
median value of the logarithm of the frequency, A log v v , the 
median of the peak duration, AT P , and the median of the 
peak frequency, v v . We also choose the highest of the peaks 
in each sequence and give for such peak its duration, AT^, 
frequency u^, and uncertainty of its localization, Alog^. 

We see that most of the peaks are practically unresolved 
since their duration and the uncertainty of the frequency is 
comparable to the Heisenberg limits of the wavelet analysis. 
For the Morlet wavelet with 2-kvq = 6 the minimum value 
of A log;/ = 0.11, and the minimum value of the time reso- 
lution is 0.24 s at 3 Hz (0.71 s at 1 Hz and 0.07 s at 10 Hz). 
The count rate of the source is high (~ 6000 cts s _1 ) which 
allows to reach the formal limit of the adopted approach. 

Searching the whole lightcurve we have found several 
interesting sequences. One of them is shown in Fig.[7| There 
seems to be a chain of peaks, systematically moving towards 
higher frequencies. Such a development could be consistent 
with propagation of the perturbations towards the gravity 
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Figure 5. (Left panels) Wavelet power spectra of Cyg X-l during its failed state transition on 1999.12.05 mapping a time-frequency 
evolution of ~ 3 Hz quasi-periodic feature (calculated for randomly selected 60 s time interval from the PCA time-series). Top values 
of wavelet power are denoted by gradual brightening of the color. Black solid contours denote significance level of 90% (outer) and 99% 
(inner) for detected peaks. 3-D visualization of the 1090—1110 s wavelet map is shown in Fig. EI A cone of influence is marked as dashed 
region. (Right panels) Corresponding global wavelet power spectra, drawn as integrated over time wavelet power spectrum multiplied by 
the frequency (see Fig. 4 for a comparison with the Fourier PSD of the same time sequences). Dotted lines represent the 95% significance 
level calculated according to the procedure given in TC98. 



center, with the expected time-scales decreasing with the 
radius. We discuss it in Section ITxHl 

For the purpose of a comparison, the contours corre- 
sponding to 99% significance level were drawn additionally 
on the wavelet maps. We found that at this level of signif- 
icance only the most prominent features of the local time- 
frequency variability are still detected. They cover a rela- 
tively small fraction of time in comparison to 90 % con- 
tours, with larger differences between the consecutive time 
sequences. Therefore, it seems that 90% contours are more 
convenient for characterizing the activity level, k, of the 
lightcurve, whereas 99% contours should be used in veri- 
fication of any trends in the frequency domain. 



5 DISCUSSION 

The nature of the short time-scale variability in X-ray emis- 
sion of accreting binaries is still under discussion. The PSD is 
generally broad-band so the variability is successfully mod- 
eled as a shot noise (Terrell 1972; Lehto 1989). Later, more 
specific models with some physical background were devel- 
oped (disk turbulence, Nowak & Wagoner 1995; magnetic 



flares in the disk corona, Poutanen & Fabian 1999, Zycki 
2002; propagation of perturbations in the inner hot flow, 
Bottcher & Liang 1999, Zycki 2003; perturbations in the 
accretion rate in the cold disk, Mineshige, Ouchi & Nishi- 
mori 1994, Lyubarskii 1997). Differentiation between those 
models is very difficult. 

Occasionally, the power spectrum shows a hint of a sin- 
gle oscillation which dominates the variability. Such a sit- 
uation happens during a transition state. Careful study of 
this apparently simple state should reveal the nature of vari- 
ations more easily. For example, during a failed transition 
state, the PSD of Cyg X-l was quite well represented by a 
single Lorentzian peak (P03). 

Such a single Lorentzian peak represents a damped co- 
herent oscillation with a time-independent frequency. This 
damping is responsible for the width of the peak. However, 
similarly broad power spectrum can in principle result from 
oscillations which are not so strongly damped but which 
change the frequency with time, like in case of a chirp signal 
(see e.g. Cohen 1995). Such a change in oscillation frequency 
may be, for example, connected with the accretion process: 
if the wave or hot plasma moves inward, the local charac- 
teristic scales may decrease with radius. 
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Figure 6. 3-D presentation of the wavelet power spectrum: (upper) for 20 s long PCA lightcurve; (bottom) for 20 s segment from a 
simulated lightcurve containing randomly occurring flares characterized by the Lorentz profile (see Section l5.H . In order of visualization 
of wavelet power distribution up to 100 Hz, maps were calculated for the lightcurves with the bin time of At = 2 -9 s. 
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Figure 7. (Central panel) Wavelet power spectrum of Cyg X-l revealing a puzzling chain of peaks ranging between 383—391 s of the 
analyzed PCA lightcurve (top panel). In the right, corresponding global wavelet power spectrum. See Section l5.8l for discussion. 
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The Fourier analysis cannot differentiate between the 
damped oscillation and the modulated frequency. Therefore, 
in order to search for the underlying physical process we 
need more advanced method of the analysis. 

In this paper we applied wavelet analysis to the Cyg X-l 
lightcurve. The wavelet maps show that strong, well local- 
ized oscillations are present rather occasionally, at frequen- 
cies ~ 3-5 Hz and they typically last for ~ 0.3 — 0.5 s. A 
few such strong peaks are always present in each of the 20 
s sequences and they effectively determine the character of 
variability, although some peaks above the 90% confidence 
threshold are typically seen for more than half of the ob- 
servation time. We did not find any systematic frequency 
change in oscillations. The analysis gave an upper limit on 
the rate of the frequency change during a single oscillation 
5lnv /d]nt < 0.2. 

This result supports the option of a coherent oscilla- 
tion with strong damping present in the system. However, 
in order to analyze more qualitatively the description of the 
system as a damped oscillator we performed numerical sim- 
ulations. 




" 0.1 1 10 



Frequency, v (Hz) 

Figure 8. Comparison of the fits to the PCA power spectrum 
performed by Pottschmidt et al. (2003) (dotted line) and provided 
by the application of equation 1101 (solid line). Diamond markers 
denote power spectrum calculated from 1000 s long simulated 
lightcurve (see Section l5.ll . 



5.1 Lightcurves simulated in time domain 

In order to model the observed lightcurve dominated by the 
Lorentzian peak we basically follow the approach of Zycki 
(2003). We create a simulated lightcurve as a superposition 
of the random shots with the shape: 

F(t) = A l£ ' (t - U)/T 



i[2nv2(t — U) + 4>i] for t > U 



(9) 



and F(t) = for t < U. The values of the frequency v 2 and 
of the characteristic time-scale r were fixed, i.e. the same 
for all shots. The amplitude of a shot, Ai, was assumed to 
cover uniformly the range between and Aq, the shot phase, 
tf>i, was assumed to be distributed uniformly between and 
27r. The moments of shot appearance, ti, were chosen from 
the Poisson distribution, for an assumed mean frequency 
of shot generation, A. It means that the time separation 
between the consecutive flares, St = U+i — ti, was given 
by 5t — — log(rand) / X where rand is a random number 
between and 1. 

It is easy to show that such a lightcurve is characterized 
by a power density spectrum: 



T(v) 



C 



+ An 2 (is — 1/2)' 



+ 



C 



+ ±n 2 (v + v 2 )' 2 



(10) 



The first term has the Lorentzian shape, as requested. 
The second term is a smooth and slowly varying func- 
tion so the sum of these two terms roughly preserves the 
Lorentzian shape. We cannot have a better representation 
of a Lorentzian power spectrum since the power spectrum 
of a real function in the time domain must have a power 
spectrum symmetric in the frequency. The Lorentzian shape 
itself is not symmetric so it corresponds to unphysical com- 
plex function in the time domain. 

This complication means that formally we cannot use 
the published values of the Lorentzian peak's parameters in 
order to generate the simulated lightcurves with the same 
power spectrum. Instead, we should refit the original power 
spectrum with the new function given by equation 1 1 Oil . 
Since this shape is not strongly different from a Lorentzian 
(see Fig. |8J, it represents the data equally well. The new 



values of the constants involved are: C = 0.945 ± 0.065, 
v 2 = 1.817± 0.184 Hz, t = 0.067± 0.003 s. New value of the 
frequency and the damping time-scale are only slightly dif- 
ferent from the values obtained by P03 from fitting a single 
Lorentzian peak (1.769 Hz and 0.07 s, correspondingly). 

The lightcurve given by equation © has statistically 
zero mean so we must add the value representing the mean 
count rate, x, in the PCA observation. The coefficients 
Aq and A must be chosen in such way that the resulting 
lightcurve has the required rms, the same as PCA data. 
Interestingly, we found that there exists a whole family of 
parameters (An, A) which reproduce the required normaliza- 
tion of the power spectrum since the constant C in equation 
HOI is given by: 



Coc 



what implies that the following relation: 
AgA = const 



(11) 



(12) 



holds for every couple of (An, A) at assumed mean count rate, 
x, and fractional rms variability. 

We also included the white noise in our simulated 
lightcurve. For this purpose, in each of the time bins of the 
simulated time-series we calculate the number of photons, 
we treat this value as 'average', draw an actual number of 
counts from the Poisson distribution around this value and 
finally we divide this value by the bin size in order to return 
to count per second units. It can be shortly written as: 



F N (t) = poiSDEv[F(i)Ai]/At. 



(13) 



Above, F(t) is given by equation @, poisdev represents a 
subroutine which generates a random number with the Pois- 
son distribution (Press et al. 1992) and thus Fjv(t) denotes 
simulated lightcurve with the Poisson noise level included. 

The level of noise does not depend on a specific choice 
of a pair (Aq, A) if it satisfies the relation 11211 . However, 
created lightcurves differ visually: low As naturally result in 
rare steep peaks while high As give apparently very noisy 
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Figure 9. (Left panels) Wavelet power spectra of the simulated lightcurve (A = 20, Aq = 2500 cts s — 1 ) calculated for randomly selected 
60-second time interval. Top values of the wavelet power are denoted by gradual brightening of the color. Black solid contours denote 
significance level of 90% (outer) and 99% (inner) for detected peaks. A cone of influence is marked as dashed region. (Right panels) 
Corresponding global wavelet power spectra, drawn as integrated over time wavelet power spectrum multiplied by the frequency. Dotted 
lines represent the 95% significance level. 



curve without clear pattern. Visually, lightcurves created 
with A ~ 20—30 look most similar to the observed lightcurve 
of Cyg X-l. We will return to this point later. 



5.2 Wavelet maps for simulated lightcurves 

We generated 1000 s long lightcurve, Fjsr(t), with assumed 
parameters: r = 0.067 s, A = 20, Ao = 2500 cts s" 1 and 
U2 = 1.817 Hz. The mean countrate of our simulated time- 
series yielded to be equal 6081.8 cts s" 1 whereas fractional 
rms 21.9%. In the first step we computed PSD to check 
whether such composed lightcurve is able to reproduce ob- 
served shape of power spectrum from the PCA data. Result 
of it has been presented in Fig. |H]by the diamond markers. 
As one can see, the shape and the width of PSD of sim- 
ulated time-series agrees very well with the main spectral 
component of Cyg X-l PSD. 

Next, from Fjv(t) we again chose randomly three se- 
quences, 20 s each, for the wavelet analysis. The wavelet 
maps for the simulated lightcurve are shown in Fig. [5] 3-D 
picture of exemplary series is shown in the bottom panel of 
Fig.il 

The maps are very similar to those obtained from the 
data analysis. This can be even better seen from the quan- 



titative parameters given in Tabled Peaks in the maps for 
simulated lightcurves also are distributed in a broad fre- 
quency range. 

At first glance this result may seem puzzling. Since 
the basic frequency in the simulated lightcurve is fixed, we 
might expect better localization of the peaks in frequency. 
However, the random choice of the phase, together with ex- 
tremely strong damping, creates an apparent diversity of the 
shapes among the consecutive shots, leading to relatively 
broad peaks, at broad frequency range, as in the data. 

Closer inspection suggests also certain systematic dif- 
ferences. Simulated lightcurves create slightly fewer peaks 
in the maps, and the median frequency in the data is some- 
what higher. This is most probably related to an excess of 
the power at high frequencies above a single Lorentzian peak 
seen in the data (~ 8-12 Hz, see Fig. 3). Otherwise, all prop- 
erties of the observed lightcurve are well represented. 

We also analyzed simulated lightcurves assuming lower 
and higher values of A (with appropriate Ao given by relation 
1121 1. Maps for significantly larger or significantly smaller 
values of A were different from the maps for Cyg X-l. Most 
noticeably, the activity level shortened to 45% for A = 5 and 
55% for A = 55; it was 59.5% for A = 20 and 61.9% in the 
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Figure 10. Wavelet power spectra of Cyg X-l corresponding to the PCA data sequence between 1090-1110 s, calculated applying 
(top panel) Morlet wavelet with 2ttvq = 3, (middle panel) Paul wavelet with m = 6 and (bottom panel) Marr wavelet. In the right, 
corresponding global wavelet power spectra. 



data maps. Therefore, wavelet maps give constraints for A, 
although they are not very sensitive to its choice. 

We tested the map sensitivity to the adopted bin size in 
the data and in the simulations. Adopting At — 2 -9 s (i.e. 
the minimal bin size available in analyzed PCA data set) 
the noise in the lightcurves (and in the power spectrum) 
increased but the maps were barely changed, as already no- 
ticed by Liszka, Pacholczyk & Stoeger (2000b). 



5.3 Maps for lightcurves simulated in the 
frequency domain 

An additional effort was undertaken in order to compare 
our simulations to the results coming from the Timmer & 
Konig (1997) approach (TK). Namely, their method allows 
to obtain a simulated lightcurve from the PSD of a given 
shape. 

Since a typical PSD of Cyg X-l is often well described 
by a power-law with two breaks and the three slopes of 0, 1 
and 2, we decided to generate first a TK lightcurve for such 
a power spectrum. The break frequencies were set at 0.1 Hz 
and 10 Hz, correspondingly, and again, randomly selected 
60 s sequence from 1000 s long time-series (At — 0.025 s) 
was selected for a closer investigation in the 1-10 Hz range. 



In result we found that the wavelet maps reveal almost no 
patterns of variability at the level of significance greater than 
90%. We estimated the n level for 20 s time-series at only a 
few per cent. Only a few peaks (1-3) were able to reach the 
significance of 99%. 

The same method and approach were used by us to gen- 
erate the power spectrum assuming the shape of the underly- 
ing variability due to a Lorentzian damping. We followed the 
numbers given in Section f5.1l to perform simulations. As an 
outcome, surprisingly we found a very good agreement of the 
wavelet peak properties comparing to the results obtained 
from the wavelet maps for the @ lightcurve. That makes 
these two techniques essentially equivalent in the sense of 
the wavelet map results. 

In none of the performed TK simulations we notice any 
chains of features with a systematically increasing frequency 
along the time axis. 



5.4 Morlet wavelet with non-standard parameters 
and other wavelet shapes 

Following the most recent analysis of De Moortel, Mundat 
& Hood (2004) we tried to increase the time resolution at 
the expense of frequency resolution by adopting smaller val- 
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ues of We constructed the Morlet wavelet map both 

for the data (see Fig. 1101 top panel) and for the simulated 
lightcurve (map not shown) using 2ttvq = 3. However, in this 
case we loose the resolution in frequency. Features become 
more elongated. Again, no pattern in frequency evolution 
can be seen. The extension of the peaks in the time domain 
decreased. A systematic difference between the observed and 
the simulated lightcurve with respect to the number of lo- 
calized peaks and the extension of active periods remains. 
Conversely, an attempt in adopting the Morlet parameter 
2nvo > 6 returned in more precise frequency resolution of 
existing features at the expanse of their time localization 
(map not shown). Also in this case, the statistics for the 
map did not change for the advantage of adopting Morlet 
function with 2-kvq > 6 2 . 

We also tried to apply two other wavelet shapes: a Marr 
and a Paul wavelet. Marr gave a very poor resolution in 
frequency both for the data fFig. 1101 bottom panel) and for 
the simulated lightcurve (map not shown), even worse than 
Morlet wavelet with 2nvo = 3. Most of the peaks were so 
elongated that they extended beyond the studied frequency 
range, making this wavelet rather useless for the studied 
lightcurves. However, it indicated a fast disappearance of the 
peaks, in time-scales shorter than 0.3 s. Also the systematic 
difference between the observed and the simulated lightcurve 
was again seen: activity periods in Cyg X-l covered larger 
fraction of time than in the simulated time-series. 

The Paul wavelet (Fig. 1101 middle panel) had effec- 
tive properties similar to the Morlet wavelet with 2-kvq = 
3. Therefore, the use of the canonical Morlet shape with 
2-kvq = 6 seems to be indeed the most profitable in the con- 
text of the X-ray lightcurves. 



5.5 Applicability of a single Lorentzian model to 
the failed transition state 

Our wavelet analysis generally supports the interpretation 
of Cyg X-l variability during the failed state transition as 
damped oscillations with fixed underlying frequency. Addi- 
tional support comes from the recent paper by Feng, Zhang 
& Li (2004) . They apply a different technique to analyze the 
lightcurve of Cyg X-l (w spectral analysis) which is sensitive 
to the shortest time-scales present in the system. The result- 
ing time-scales for the hard and the intermediate state are of 
the same order as the damping time-scale of the Lorentzian, 
consistent with the extension of the detected features on our 
wavelet maps. 

Simulated lightcurves gave maps quantitatively similar 
to the maps obtained from the Cyg X-l lightcurve. The ob- 
served PSD contains an excess of power at higher frequencies 
(modeled by P03 as an additional Lorentzian peak). Simu- 
lated lightcurves did not contain this additional high fre- 
quency input but still reproduced the maps properly so the 
dominant Lorentzian is responsible for the character of the 
variability both in the frequency domain and in the time- 
frequency domain. 



2 Please note that Morlet wavelet with 2itvq < 5 and > 6 does 
not meet admissibility condition ('Scction l3.2.1l and its usage can 
be justified only for non-quantitative analysis. 



5.6 Dominating PSD components in hard, soft 
and failed soft state 

Our results give a strong support to the interpretation of 
the PSD through Lorentzian peak during the failed state 
transition. It means that in general case we most probably 
deal with three main physical components of the PSD with 
the following properties (Nowak 2000; P03): 

• Power law component v~ : this component dominates 
the PSD at v < 10~ 5 Hz (Reig, Papadakis & Kylafis 2002), 
it extends up to 15-20 Hz when the source is in soft state, it 
extends into the region v ~ 10 _3 -10 Hz but with suppressed 
amplitude during the transition state and it practically dis- 
appears in the hard state; 

• Lorentzian peak L%: present in the hard state but its 
amplitude strongly decreases during the transition state, ac- 
companied with the rise of the power law in high frequency 
range; 

• Lorentzian peak L2: suppressed only when the the soft 
state is reached; ratio of the two Lorentzian frequencies is 
constant (~ 7) although the frequency can change in long 
time-scale up to a factor of 4. 

Results for the coherence function and the Fourier time 
delays support the view that both Lorentzians are generally 
highly coherent, also between themselves (Nowak et al. 1999) 
but when one of the peaks disappears in the failed state 
transition, the mean time delay between 2-4 keV and 8- 
13 keV increases to 8 ms (a typical value if measured at ~ 
3 Hz), and the coherence in 3-10 Hz band drops to 0.6 (P03). 
Similar properties were observed by Cui et al. (1997b) (e.g. 
30.05.1996) when the source was on its way to the completed 
transition to soft state. 

These components seem to be critical for our under- 
standing of the nature of the observed variability although 
exact fits to the data generally require more components. 



5.7 Physical interpretation of the X-ray 
variability 

5.7.1 Power Law component v~ x : propagating disk 
oscillations filtered at the transition radius 

There are several independent arguments in favor of the hy- 
pothesis that the v~ x power law component of the PSD is 
related to the presence of the cold disk. This component 
is exclusively responsible for the observed variability when 
Cyg X-l is in the soft state, i.e. radiation spectrum is dom- 
inated by the disk emission and we see the relativistically 
broadened iron Ka line which gives constraints for the cold 
disk inner radius (e.g. Di Salvo et al. 2001). The time-scales 
involved cover uniformly a large range, including time-scales 
as long as years which are naturally expected in the outer 
disk. However, there were two puzzles which were not ini- 
tially understood. First, what we see is not the variable disk 
emission, since the direct disk emission is roughly unchanged 
but variable Compton component (Churazov, Gilfanov & 
Revnivtsev 2001; Maccarone & Coppi 2002). Second, the 
energy contained in the longest time-scale variability is con- 
siderable while the gravitational energy at the outer edge of 
the disk is small. 

The model which successfully explains the situation is 
the scenario of propagating perturbations. The idea was 
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nicely outlined by Lyubarskii (1997), although some basic 
ideas were already contained in the automaton model of Mi- 
neshige et al. (1994). More sophisticated approach was devel- 
oped by King et al. (2004). In this model local perturbations 
form in the disk as a result of the action of the magnetoro- 
tational instability (MRI; Balbus & Hawley 1991). Charac- 
teristic frequencies of these perturbations at a give radius r 
are somewhat lower than the local Keplerian frequency and 
can be expressed as: 

--H£)(^rm » <»> 

Most of the authors argued that the inner radius of the disk 
in the soft state is of order of 3 -Rschw (e.g. Frontera et al. 
2001 give fttgRg, where R g = 0.5i?s c h w ) and we adopted 
this value as a convenient unit. The characteristic frequency 
at this radius reproduces the observed PSD frequency break 
at 20 Hz in the soft state when fc mag = 5.5. The factor of 
2n was included, i.e. v — Q/(2tt). Those perturbations prop- 
agate inward, to the region where dissipation and the pro- 
duction of the hard X-ray emission proceeds. Therefore, the 
accretion rate at the inner dissipation region has a memory 
of all time-scales present in the disk at all radii. 

When the source proceeds from a soft state to a hard 
state, the power law component is not seen at time-scales 
shorter than 10 _5 -10 -3 Hz. We cannot relate this change 
to a change of the disk inner radius since this would re- 
quire a change of the radius by a factor of 10 3 -10 4 . Analysis 
of the data suggest much lower values; e.g. ~ 5-35i?s c hw (di 
Salvo et al. 2001), 10±|i^ (Frontera et al. 2001) or ~ 100i? g 
(Gilfanov, Churazov & Revnivtsev 2000). However, due to 
a different treatment of the Compton thickness of the inner 
disk, the values of the radii quoted above should be consid- 
ered here with a caution. 

However, a natural explanation of this change is found 
if we adopt the strong ADAF principle after Narayan & Yi 
(1995). This is a statement which tells that whenever an 
ADAF can form, it does form. The consequence of this as- 
sumption is that for a given value of the transition radius, 
r tr , only one value of the accretion rate is possible. This fact 
automatically leads to filtering all time-scales in the momen- 
tary accretion rate shorter than a local viscous time-scale. If 
the perturbed accretion rate is somewhat larger than allowed 
by a current position of the transition radius, the matter ac- 
cumulates near r tr while the accretion rate in the inner flow 
is unperturbed. Mass accumulation, if persisting for long 
enough, leads to a build-up of a cold disk in part of the re- 
gion previously occupied by the inner hot flow. The position 
of the transition radius moves in, and after this change an in- 
ner hot flow can transmit mass at higher rate, appropriately 
to the new position of r tr . Alternatively, if the perturbed ac- 
cretion rate in the disk is somewhat lower than current rt r , 
the inner hot flow still persists with the same accretion rate 
at the expense of the cold disk. Finally, the cold material is 
used up, the transition radius moves out to the new position, 
and the accretion rate in the inner flow decreases in agree- 
ment with new r tr . Since both the removal and the buildup 
of the cold disk happens in a local viscous timescale all vari- 
ations in the accretion rate slower than that are transmitted 
while faster variations are damped. The fastest timescale for 
a transition from a hard to a soft state is also given by a 



viscous time-scale, 

— L5xl ° (— ){—) (so^J x 

and it is observationally constrained to be of order of a day, 
or a fraction of a day, consistent with the time-scale 10 -3 - 
10" 5 Hz. Here we assumed the expected ratio of the disk 
thickness to the radius, h/r of order of 0.08 and a = 0.01 
after King et al. (2004), for the scaling purposes. 

This picture gives a strong observational support to 
the strong ADAF principle. An independent argument for 
this kind of behavior was also found by Czerny, Rozariska & 
Kuraszkiewicz (2004) from the study of the constraints for 
an inner radius from Broad Line Region in active galactic 
nuclei. 

During the transition state the accretion rate is most 
probably much more strongly enhanced than during the 
hard state perturbations. Observations show a sudden de- 
crease in the inner radius. For example, on 1996.05.30, when 
the state was yet similar to the failed state transition, the 
inner radius was estimated to be 9l4 3 i?schw (Gierliriski et al. 
1999). However, we can suspect that the material piles up 
rapidly. Instead of smooth motion of the transition radius 
clumps of the cold material possibly enter the inner hot flow 
region and accrete. Such a direct accretion of a cold phase 
may (i) explain why during the transition, including failed 
transition state, the short time-scale perturbations in the 
accretion flow start to propagate inward, (ii) cold moving 
blobs provide additional source of soft photons for Comp- 
tonization which disrupt the high coherence seen during the 
hard state. 

It still remains to explain what kind of instability leads 
to the disk destruction when ADAF solution is possible. An 
interesting scenario is discussed by Spruit & Deufel (2002), 
but it requires the initial existence of the hot ions. It may 
give a hint why we observe a strong hysteresis effect in 
some sources but not in Cyg X-l. This hysteresis (see e.g. 
Miyamoto et al. 1995; Meyer-Hofmesiter, Liu & Meyer 2005; 
Zdziarski et al. 2004) leads to a transition from soft to hard 
state at much lower accretion rate than the transition from 
hard to soft state. In the soft state of Cyg X-l significant 
fraction of emission is of non-thermal origin. Hot ions are 
always there, and the mechanism of Spruit & Derfel (2002) 
may operate while we need another mechanism for sources 
which do not show any hot plasma in the soft state. 

5.7.2 Lorentzian peaks L\ and L2: MRI and dynamical 
pulsations of the inner ion torus 

The two Lorentzian peaks dominating the hard state are 
expected to be related to the inner ion torus. The constant 
ratio between their frequencies strongly suggest that we see 
two kinds of variability coming from the same medium. 

Since the inner hot flow must transport the material the 
same MRI instability is expected to be operating. Since the 
inner flow is expected to be geometrically thick, instability 
has more global character there, and the dominant frequency 
is likely to be determined by the transition radius. This fre- 
quency still should be roughly reproduced by equation 1141 : 
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When the second Lorentzian peak is at 1.7 Hz, the first 
peak is at 0.25 Hz (P03), and such a frequency is repro- 
duced for a reasonable value r tr = 55i?schw In general, the 
position of the Lorentzian varies. Revnivtsev, Gilfanov & 
Churazov (2001) showed that the corresponding peak in the 
power spectrum of GX 339-4, described as quasi-periodic os- 
cillation (QPO), strongly correlates with the other spectral 
properties of this source, supporting its connection with the 
change of the transition radius. 

Higher frequency variations are likely to represent the 
dynamical pulsation of the inner hot flow, or traveling sound 
waves. Global oscillations of this medium, modeled as a 
torus, were studied in a number of papers (e.g. Abramow- 
icz, Calvani & Nobili 1983; Giannios & Spruit 2004; Mon- 
tero, Rezzolla & Yoshida 2004). The typical frequencies of 
p-modes depend on the details of the model (e.g. assumed 
radial distribution of the angular momentum) but they are 
roughly of the order of the Keplerian frequency at the outer 
edge of the torus (Rezzolla, Yoshida & Zanotti 2003). The 
excitation of the modes is due to a radiative interaction of 
the disk and the torus, while the synchrotron emission pro- 
vides the damping mechanism (Giannios & Spruit 2004). 

The expected frequency of this mode is roughly: 
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so the value of the transition radius ~ 50i?schw gives the 
value of the frequency roughly corresponding to the observed 
frequency of the Lorentzian, 1.7 Hz. 

During the failed transition state the first Lorentzian 
peak is strongly damped, and a power law components is 
partially rebuilt, as discussed in Section 15.7.11 If the first 
Lorentzian peak is connected with MRI instability, the ex- 
change of power between a power law component and the 
first Lorentzian is natural. Decreased MRI instability in the 
inner flow suppress the accretion of the hot material and 
the accretion proceeds predominantly through cold blobs. It 
is an interesting question whether the presence of the cold 
blobs prevent MRI instability or the absence of hot material 
accretion forces the cold phase to take over but the even- 
tual compensation of the two mechanisms is quite intuitive. 
Apparently, dynamical pulsations persist as long as the hot 
inner torus is still there. If a transition to a soft state is 
completed, this phase finally disappears. 

The suggested picture may be too simple since actu- 
ally a number of instabilities may operate in the accretion 
flow (e.g. Menou, Balbus & Spruit 2004 and the references 
therein) but we consider our interpretation as a plausible 
and attractive scenario. Clearly, quantitative models would 
be needed to work out the detailed predictions. 

Observationally, an insight may come from the pres- 
ence of rapid, very energetic events reported by Gierlihski 
& Zdziarski (2003) both during soft and hard state. It is 
interesting that a soft state flare no. 13, studied in detail 
by these authors, was found to be well described by damp- 
ing time-scale (after the peak) of r = (21 ± 6)°- 7±0 1 m s 5 
i.e. t ~ 0.07 s. This damping time-scale is the same as the 
damping time-scale of the L2 Lorentzian, which may indi- 



cate that damping is related to magnetic field reconnections. 
Same event observed in the hard state yielded much longer 
damping time-scale, ~ 0.32 s which in turn roughly coincides 
with the damping time-scale of the first Lorentzian, L\. Al- 
though the exact origin of the flares is not known, one takes 
into account sudden conversion of energy accumulated in the 
disk into magnetic heating of a hot plasma or fast release of 
energy due to magnetic field reconnection in the flares hung 
above the disk (Gierlihski & Zdziarski 2003 and references 
therein). The difference between the hard state and the soft 
or failed transition state may be due to the absence of the 
cold material in the former case. Therefore, in the hard state 
the large scale magnetic loops exist within the torus, while 
in the soft state or failed transition state magnetic field lines 
are partially frozen into the cold disk or cold blobs. 



5.8 Wavelet analysis as a tool for tracing 
accreting matter onto black-holes? 

It is tempting to associate the observed puzzling chain of 
peaks in Fig. |7| (mentioned already by us in Section 14. H 
with a systematic trend. Closer look in this figure shows up, 
in fact, not one but two sequences of peaks roughly of the 
same slope. The first one, most prominent, extends between 
383-391 s, whereas the second one between 389-397 s. Com- 
parison of wavelet map with corresponding PCA lightcurve 
reveals that all wavelet features are associated with some 
peaks occurring in the same moments in the time-series, 
as expected. In addition, a careful inspection of all 20 s 
long sequences of the PCA observation uncovers more simi- 
lar chains of the same slope. Unfortunately, their number is 
low i.e. ~ 14 per whole lightcurve, each one of duration < 7 
s. Interestingly, we did not find similar chains of evolution 
from lower to higher frequencies in the simulated data. 

The time-frequency evolution of some features is not to- 
tally unexpected. Both the hot ion torus material and cold 
clumps flow onto the black-hole, in the direction of decreas- 
ing Keplerian time-scales. We can witness a manifestation 
of the single events like a sequence of magnetic field recon- 
nections in an inflowing material. Theoretically, this issue 
were discussed already by Stoeger (1980), however, his con- 
siderations were mostly appropriate for events taking place 
below the marginally stable orbit. Here, it seems that ob- 
served features between 1-10 Hz correspond rather to the 
range of 69-15 -Rschw • 

If these enigmatic chains of oscillations are related 
somehow to the processes of X-ray emission coming from 
the material accreting onto black-holes, then, in our opin- 
ion, a wavelet analysis may become a promising tool in their 
spotting. However, any judgement on the reality and phys- 
ical origin of the observed trends can be only done after 
much more quantitative analysis. The number of the de- 
tected chains must be large and statistically significant. 



6 CONCLUSIONS 

For the first time, we applied the wavelet multi-resolution 
approach in order to perform the complex studies of aperi- 
odic X-ray variability of Cyg X-l. Considered data set cor- 
responded to its failed state transition on 1999.12.05. On 
that day, the Fourier power spectrum revealed a broad-band 
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peak between ~1-10 Hz, well described by a single Lorentz 
function of unclear origin. On the contrary to the Fourier 
method, the wavelet analysis provided an excellent inspec- 
tion of the frequency content of the signal as well as allowed 
for simultaneous spotting the time evolution of the extracted 
features. 

Our results can be summarized as follows: 

• Cyg X-l lightcurve analyzed in the time-frequency 
plane contained randomly appearing flare-like structures of 
duration ~ 0.5 s and of the average median of the peak 
frequency width Alogt'p ~ 0.14; 

• observed Lorentzian shape forms through summation of 
these time-localized structures; 

• most of the Cyg X-l lightcurve properties are well re- 
produced by the simulated lightcurve composed of a series 
of randomly occurring flares. In our model, each flare was 
described as a damped oscillator with variable amplitude 
and phase but fixed excitation frequency and time-scale of 
decay; 

• Cyg X-l wavelet maps may show insignia of an inflow 
while the maps representing the simulated lightcurve do not 
seem to show similar tendency. However, no statistical sig- 
nificance was assigned to this hypothesis. 
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APPENDIX A: DISCRETE FORM OF CWT IN 
TERMS OF FOURIER TRANSFORM 

Let us express a continuous wavelet transform as the convo- 
lution of the Fourier transforms of the lightcurve and wavelet 
function in the following way: 



Wa,b = 



The Fourier transform of ip a ,b is given by: 

j> a , b (v) = 4= r V- (—} e-^dt. 
\/a ./ ^ V a 



(Al) 



(A2) 



Substituting (t — b)/a by r\ equation IA2I evaluates into: 

r oo 

V'a.bM = Vae~ l27T " b / ii(r/)e~ t27ra ' /r 'dri (A3) 



where a constant component was separated out. As one can 
see, an integral expression can be considered as the Fourier 
transform of a wavelet function where 2nav is a rescaled 
angular frequency. Thus it is possible to denote IA3I shortly: 



and therefore a CWT as the inverse Fourier transform of the 
product: 



W a ,b 



x(v)xj)* {2irai>)e 2nl ' 1 dv 



(A5) 



where complex conjugate of the Fourier transform of tp func- 
tion was taken into account. 

A discrete Fourier transform of x k is given by 



-i2irjk/N 



(A6) 



where j denotes a frequency index. Performing a discrete 
notation of the time-shifting factor as b = kAt where At 
stands for the sampling time of a lightcurve, one can write 
down a discrete form of continuous wavelet transform in 
terms of the Fourier transform as: 



w k (a) 



f2na\ 
VAT) 



1/2 - 



Xj^*{2-Kav ] )e l 



3=1 



where frequency Uj equals: 

j/(N ohB At) for j < N ohB /2 
-j/ (N oba At) for j > N obs /2 



(A7) 



(A8) 



and the factor (2na/ At) 1 ^ 2 ensures that wavelet will keep 
the same energy at every scale a. Because of the finite dura- 
tion of the signal, the proper choice of scales must be done. 
The largest scale would correspond to the lightcurve span 
of T whereas the smallest one to ought to be an equivalent 
of the Fourier Nyquist frequency, i.e. ao = 2At. In fact, a 
selection of scales in this region can be chosen arbitrarily 
with a step in scale not smaller than ao. However, in case of 
very fine signal when a few decades of scales are wished to 
be covered, one can build the wavelet transform choosing: 

-»mAm 



0,...,M 



(A9) 



dm — a 0<Z 
where 

M = Am" 1 log 2 (N ohs Ata,o 1 ) (A10) 
(TC98) and thus rewriting <A7> in the following way: 



uik(am) 



2Tva„ 



V At 



1/2 JVobs 



Xjip* (2-Ka m Vj)e l 



(All) 



This paper has been typeset from a TpX/ ETpX file prepared 
by the author. 
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